Non-local modulation of the energy cascade in broad-band forced turbulence 
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Classically, large-scale forced turbulence is characterized by a transfer of energy from large to 
small scales via nonlinear interactions. We have investigated the changes in this energy transfer 
process in broad-band forced turbulence where an additional perturbation of flow at smaller scales 
is introduced. The modulation of the energy dynamics via the introduction of forcing at smaller 
scales occurs not only in the forced region but also in a broad range of length-scales outside the 
forced bands due to non-local triad interactions. Broad-band forcing changes the energy distribution 
and energy transfer function in a characteristic manner leading to a significant modulation of the 
turbulence. We studied the changes in this transfer of energy when changing the strength and 
location of the small-scale forcing support. The energy content in the larger scales was observed to 
decrease, while the energy transport power for scales in between the large and small scale forcing 
regions was enhanced. This was investigated further in terms of the detailed transfer function 
between the triad contributions and observing the long-time statistics of the flow. The energy is 
transferred toward smaller scales not only by wavenumbers of similar size as in the case of large-scale 
forced turbulence, but by a much wider extent of scales that can be externally controlled. 

PACS numbers: 47.27.E-, 47.27.Gs, 47.27.Rc 



I. INTRODUCTION 



The dynamics of kinetic energy plays a central role 
in turbulent flows. The nonlinear term in the Navier- 
Stokes equations is responsible for the transfer of en- 
ergy between any three wavevectors that form a triad 
in spectral space [l|. Along with the viscous and forc- 
ing terms this controls the production, transfer and dis- 
sipation of energy in the system. The triadic interac- 
tions have been studied for decaying and forced turbu- 
lence by many authors (for a review see |2|). Through- 
out the years various types of large-scale forcing meth- 
ods HHSE0EI3E3II3 have been p r °p° sed to 

sustain quasi-stationarity in numerical turbulence as an 
idealized form of turbulent flow. The aim of such nu- 
merical experiments was to investi gate the basic concept 
of the Kolmogorov (K41) theory [l^] that proposes an 
inertial range in the kinetic energy spectrum and local 
transfer of energy within this range. The turbulent ki- 
netic energy is on average transferred locally from larger 
to neighboring smaller scales. 

The purpose of this paper is to numerically investi- 



*Electronic address: a.k.kuczaj@utwente.nl 

tAlso: Anisotropic Turbulence, Fluid Dynamics Laboratory, De- 
partment of Applied Physics, P.O. Box 513, 5300 MB Eindhoven, 
the Netherlands 



gate the processes associated with the flow of energy in 
a turbulent flow. Specifically, we consider modulated tur- 
bulence in which the modifications involve the supple- 
mentary forcing in a wide range of modes located in an 
inertial range of the flow. In the literature, mainly turbu- 
lence with forcing restricted to the large scales has been 
examined in detail |2|. The small scale behavior was 
found to be energetically quite insensitive to the type of 
forcing and at sufficiently high Reynolds numbers a well- 
developed inertial range was observed [13j . Against this 
background, we extend the use of forcing methods and 
investigate their application directly in the inertial range, 
thereby focusing particularly on the competition between 
transfer and forcing. We quantify the dominant alter- 
ations due to the broad-band forcing in terms of changes 
in the energy cascading processes. We pay attention to 
the energy transfer function and consider changes that 
arise in the contributions from 'local', 'non-local' and 
'distant' triadic interactions. Compared to traditional 
large-scale forced turbulence, we observe a strengthening 
of the contributions of non-local interactions, leading to 
a modification of the inertial range spectrum. 

High-resolution direct numerical simulations of tur- 
bulence that measure the influence of individual terms 
in the Navier-Stokes equations on the triadic interac- 
tions have been reported dEHElHHIIi. It 
was found that the energetically dominant triadic inter- 
actions involve sets of three modes in which the mag- 
nitude of the wavevector of one of the modes differs 



considerably from the other two. This suggest that 
statistics of smaller scales may be affected by larger 
scales. These dominant processes are not in contradic- 
tion with the Kolmogorov theory because the energy is 
mainly exchanged between the two modes of quite similar 
wavevector-size 19]. Only a small net energy transfer to- 
ward larger wavenumbers arises that involves a detailed 
cancellation between many individual triad transfers [20] • 
The spectral space dynamics is characterized by a multi- 
tude of separate transfer-processes among various modes. 
These contributions can be collected in pairs with op- 
posite sign and almost the same magnitude. In total, 
this leads to a large number of 'near-cancellations' and 
hence only a comparably small net effect remains that 
constitutes the well-known 'downward cascading' toward 
higher wavenumbers in spectral space. This was con- 
firmed with the use of helical mode decomposition in [2(J • 

The dynamics of actual turbulent flows seen in nature 
is usually characterized by an enormous number of inter- 
acting scales, often perturbed by geometrically complex 
boundaries and influenced by additional forces such as ro- 
tation and buoyancy. This can lead to inhomogeneity and 
anisotropy, which are not covered directly in the classical 
view of the Kolmogorov energy cascade and may express 
themselves in non-local interactions of various particular 
scales of motion. The complexity of such systems mo- 
tivated us to study in more detail forcing methods that 
simultaneously perturb a prescribed range of scales [2l| . 
Such 'broad-band' agitation of various scales of motion is 
observed experimentally in turbulent drag reduction by 
fibre suspension |22.l23f , flo ws through porous media J2J] 
and over tree canopies [2j|. In these cases the energy is 
transferred abruptly to small scales when the flow reaches 
an obstruction. Various other types of flows also exhibit 
turbulent motions that coexist at different scales J2g|. 

To explore the possibilities of a broader application 
of forcing methods in turbulence modeling and concur- 
rently examine the energy dynamics in flows that do not 
directly follow the classical Kolmogorov —5/3 scaling we 
employ numerical simulations of broad-band forced tur- 
bulence. The forcing studied in this paper represents a 
continual addition/removal of energy from a broad range 
of scales in the system, thereby providing the possibility 
of altering the characteristic —5/3 slope in the kinetic 
energy spectrum as predicted by the K41 theory. Specif- 
ically, as indicated in Fig. ^ we apply the forcing to two 
regions. The large-scale forcing k < k classically agi- 
tates the largest scales in a flow while the additional band 
k\ < k < ki is located in a region of the inertial regime, 
to allow a direct competition with the nonlinear transfer 
term. For incrtial-range scales broad-band forcing intro- 
duces explicit energy injection next to the transfer-term. 
We varied the spectral support and strength of the high-fc 
band to investigate the modulation of the turbulence that 
develops. This distinguishes it from the classical forcing 
of large scales only. 

In this paper we compute changes in the energy dis- 
tribution associated with the broad-band forcing and ob- 
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FIG. 1: Broad-band forcing in spectral space. Large-scale 
forcing k < ko with an additional high-fc forced region 
fci < k < &2- 



serve a characteristic alteration in the spectral energy 
transfer compared to the classical Kolmogorov cascad- 
ing. This alteration expresses itself by additional local 
minima and maxima in the transfer function. It is well 
known that in cases with large-scale forcing only, negative 
values are found for the transfer at the smallest wavenum- 
bers indicating the energy injection at these scales. The 
positive values for the transfer that arise for all other 
wavenumbers indicates the energy cascading process to 
smaller scales. In our case of broad-band forcing in the 
inertial range, additional negative regions appear in the 
transfer function. These coincide with the additional lo- 
cal injection of energy. Such a negative region is bordered 
by nearby additional maxima in the transfer. These char- 
acterize the associated increased energy transfer to scales 
just larger or just smaller than the broad-band forced re- 
gion. 

Forcing applied to different spatial scales simultane- 
ously allows a non-local modulation of the energy distri- 
bution compared to the reference Kolmogorov case. To 
quantify the alterations in the energy transfer we use a 
decomposition of the velocity field closely following [l(| 
and investigate the magnitude of the contributions from 
various spatial scales to the overall energy transfer. 

The main finding of this study pertains to the role of 
broad-band inertial range forcing in modifying the nat- 
ural energy cascading process. This is understood ex- 
plicitly in terms of changes in the detailed non-local en- 
ergy transfer. In addition, we illustrate and quantify the 
mechanism of enhancement of the total energy transfer 
to smaller scales arising from broad-band forcing and the 
depletion of the energy-content in the large scales. Agi- 
tation of certain high wavenumbers can affect well sepa- 
rated low wavenumber components in a flow. These find- 
ings may be relevant for problems that involve the control 
of turbulent flow in complex geometries in which vari- 
ous scales of motions are simultaneously agitated, e.g., 



in compact heat-exchangers 24]. Further applications 
of such broad-band forcing may be connected with the 
observed modulations of transport properties in physical 
space leading to an enhanced scalar dispersion rate |2l| . 
The organization of this paper is as follows. The math- 
ematical formulation of the problem is given in Sec. |n] 
where the computational method and the energy transfer 
terms are also described. The energy spectra of broad- 
band forced turbulence and the modulation of the energy 
transfer are investigated in Sec. lIHI In Sec. lIVI we present 
a more detailed view of the energy transfer processes by 
computing partitioned energy transfer function over var- 
ious spatial scales. The paper closes with a summary in 
Sec. El 



II. COMPUTATIONAL FLOW MODEL 

An overview of the computational model is given in 
this section (JII Afl . The forcing method in the broad-band 
context is described subsequently (|IIB(I . In addition, the 
energy transfer between different triads is partitioned in 
a number of contributions IjII Cfl that will be studied nu- 
merically in Scc. lIVI Finally, simulation details are given 
(ITTDl). 



A. Equations of motion 

The incompressible Navier-Stokes equations in spec- 
tral (Fourier) representation can be written as 

(d t + uk 2 ) uc(k,t) = tf Q (k,t) + F Q (k,t), (l) 

where u a (k, t) is the velocity field coefficient at wavevec- 
tor k (k = |k|) and time t [l|. The non-dimensional 
kinematic viscosity v is the inverse of the computational 
Reynolds number (Re — 1/f). The nonlinear term reads 



* a (k,t) = Af a/ s 7 J2 u (P>*K(<1>*)> ( 2 ) 

p+q=k 



and the forcing term F a (k, t) is specified in section III Til 
The tensor M a/ g 7 in (J2J accounts for the pressure and 
incompressibility effects: 



in which 



M a/3l = — (kpDay + k 7 D a f3J , 



D a fj = S a fj - k a kp/k 2 



(3) 



(4) 



Taking the inner product of Q) and u*(k, i), where the 
asterisk denotes the complex conjugate, we obtain the 
energy equation 

(d t + 2vk 2 ) E(k, t) = T(k, t) + T F (k, t). (5) 



The spectral energy density is denoted by E(k, t) — 
\u* a (k, t)u a (k, t). The rate of energy exchanged at 
wavevector k with all other modes in the system is char- 
acterized by the energy transfer function 



T(k,t)=<(M)tt a (M). 

The rate of energy provided by the forcing term is 



(6) 



T F (k,i)=<(M)i^(k,t), (7) 

and the energy dissipation rate present in © reads 

e(k,t) = 2isk 2 E(k,t). (8) 

The three terms T(k, t), Tp(k, t) and e(k, t) represent 
the energy dynamics in the system that each typically 
act in distinct wavenumber regions. The forcing term 
Tp(k, t) is non-zero in the forced modes only. In this pa- 
per the collection of forced modes will always contain a 
low wavenumber band corresponding to large-scale forc- 
ing of the flow. In addition, higher wavenumber contri- 
butions will be included in Tp (k, t). In contrast, the en- 
ergy dissipation rate e(k, t) is defined in the entire spec- 
tral space, but it is dynamically important primarily for 
the high wavenumber range. Finally, the transfer term 
T(k, t) is basic to the development of an energy cascade 
and is a dominant contribution for wavenumbers in an 
inertial range [l| . 

The change of the total energy E in the system is con- 
nected with its viscous dissipation and the total effect of 
the forcing. In fact, introducing 



we find 



E(t)=T E ( k > t )> 



d t E(t) = T F (t) - e(t), 



(9) 



(10) 



where e(t) = £k e (M) and T F (t) = £k T iKk,i). We 
used the fact that the total energy transfer T(t) — 
J2w T(k, t) = 0. The injection of energy occurs only in 
the forced region. This keeps the whole system in a quasi- 
stationary state. Normally, the forced region is restricted 
to the largest scales in a flow represented by the smallest 
wavenumbers 5, 10] . The energy introduced in the large 
scales is transferred to smaller scales and dissipated pri- 
marily in very localized flow-features of viscous length- 
scales. By the introduction of an additional source of 
energy in the inertial range we will study the perturba- 
tion of the energy cascading process by the forcing. The 
forcing method adopted here will be presented next. 



B. Forcing method 

Forcing is achieved by applying an additional driving 
F a (k,t) to the velocity field in Fourier space, cf. QJ. 



Conventionally, the turbulent cascade develops as a sta- 
tistical equilibrium is reached, characterized by the bal- 
ance between the input of kinetic energy through the 
forcing and its removal through viscous dissipation. In 

literature (|lSHE0Eli[lS[I3),we may distin g uish 
several numerical approaches to forced turbulence that 
all refer to the agitation of the largest scales of motion. 
Here, we modify such classical forcing procedures by al- 
lowing for the simultaneous agitation of a broader range 
of intermediate-fc modes as depicted in Fig. ^ 

We study two ranges of forcing: the classical large- 
scale forcing (k < ko) and small-scale forcing localized in 
the spectral region where the transfer of energy T(k, t) 
is important (fci < k < fe). By narrowing or widening 
the width of the forced bands, along with a change in 
their location in spectral space we can control several as- 
pects of the energy-dynamics. The strength of forcing is 
controlled by the amount of energy introduced to various 
regions in spectral space. 

We expect the small-scale forcing band to influence 
the inter-scale energy transfer process not only between 
scales of similar size but at a wider spectrum of scales. 
This may be understood globally as follows. The process 
of energy cascading is mainly interpreted via the result- 
ing local transfer of energy in spectral space 19] . How- 
ever, this total energy transfer results from many non- 
local contributions and these may be directly altered by 
the additional small-scale forcing. Correspondingly an 
influence on the overall energy cascading process may 
occur over an extended wavenumber range. We quantify 
this effect by evaluating the nonlinear interactions among 
the various modes while they are being perturbed by the 
broad-band forcing. 

In this paper we adopt the recently proposed fractal 
forcing |27j . which involves a power-law dependence of 
F a on the wavenumber: 



The vector e in l|llf) is given by [27 



F a (k,t)=e w (k) 



fc /3 e Q (k,i) 



Eu^V2E(kJy 



(11) 



where the coefficient j3 — Df — 2 is connected with the 
fractal dimension Df of the stirrer and e w (k) is the en- 
ergy input rate at mode k. The set of forced modes IK 
is composed of bands ^ m ,p (jn < p) which consist of 
p — m + 1 adjacent spherical shells B n = §~(n — 1/2) < 

k| < 7~(n + 1/2): m < n < p, in terms of the size of 
the computational domain denoted by L\,. In the sim- 
ulations we always force the first shell Si and a single 
high-fc band K mp , if not stated otherwise. The classical 
large-scale forcing of the first shell Si has a constant en- 
ergy injection rate s w> i in l|ll(l while K mjP has a constant 
strength e w ^ and a support in spectral space controlled 
by ?n and p: 



e w (k) 



£w,l 


if keSi, 


&w,2 


if k e Km >p , 





otherwise. 



(12) 



e(k,t) 



u(k,/j) kxu(k,/j) 



+ i 



u(M)| |k||u(k,t)| 



(13) 



This vector consists of two parts, either parallel or per- 
pendicular to the vector u(k, t). In this forcing proce- 
dure, we have control over the energy input rate, the 
range of forced modes and the effective geometrical com- 
plexity of the stirrer represented by the fractal dimension. 
The summation over all forced modes of 
u* a (k, t)F a (k., t) yields a total energy input rate given 
by: 

(14) 
where e w — e w s + e w ^- The energy input leads to a 
quasi-stationary state described by the energy equation 



d t E(t) 



m- 



(15) 



This characterizes the energy dynamics in the system at 
the most global level. We observe that this forcing im- 
plies a constant energy injection rate that results in a 
fluctuating total energy E and a fluctuating total energy 
dissipation rate with mean e w . 

In the next subsection a more detailed description of 
the energy dynamics will be given based on the processes 
that are responsible for its transfer. 



C. Energy transfer 

A detailed investigation of the energy transfer in large- 
scale forced turbulence 19] shows that the dominant tri- 
adic interactions occur between wavevectors of quite dif- 
ferent lengths. Hence, large-scale forcing may be directly 
involved in the dynamics of much smaller scales [28j . The 
interactions are roughly classified as "local" when the 
sizes of all wavevectors in a triad are similar, "non- local" 
when the scale separation is about a factor 10-15 and 
"distant" when the separation is much larger [23 . It was 
shown that the transfer of energy reaches maximum val- 
ues for triads with two wavevectors of similar size and 
one with quite different length [H|. Although, the inter- 
actions between triads can be seen mainly as non-local, 
the dominant net ener gy t ransfer is local, i.e., occurring 
between similar scales [lj, LMJ, lLLI • The interactions pro- 
duce forward and backward energy transfer that com- 
bined result in a small net forward energy transfer be- 
cause of the detailed balance between contributions that 
virtually cancel each other [2(J]. The forward cascade in 
the inertial range was found to be dominated by local 
and non-local interactions, while the distant interactions 
do not significantly transfer energy |2J|. All these find- 
ings concern the classical turbulence forced at the largest 
scales. 

Against this background, we ask what the turbulence 
response will be to a broad-band perturbation of the en- 
ergy transfer processes? In recent literature a somewhat 




wavevectors p and q: 



FIG. 2: Schematic triadic interaction that occurs between 
wavevector k in shell Sfc and wavectors p, q taken from regions 
P and Q of spectral space that each consist of four shells with 



related study was reported in [33]. Decaying turbulence 
that starts from an initial condition with an energeti- 
cally strongly enhanced small-scale band of modes was 
studied. The presence of the extra small-scale band was 
found to reduce the intensity of the developing turbu- 
lence by enhancing the non-local energy cascade directly 
towards smaller scales. This removes the kinetic energy 
more efficiently. The energy feeding mechanisms and en- 
ergy transfer also attract much attention in transitional 
and turbulent flows with an active control 33] . The mod- 
ulation induced by the broad-band forcing has its conse- 
quences not only in the spectral space dynamics of a flow 
but also in its physical space transport properties |21j . 

To analyze the response of turbulence to the additional 
broad-band perturbation in more detail we apply pre- 
viously developed methods used in the examination of 
energy transfer in large-scale forced turbulence 16]. Re- 
ferring to Fig. [2 the energy transfer between a wavevec- 
tor k = (ki,k2,ks) and all pairs of wavevectors p and 
q = k p with p, q chosen in some prescribed regions 
P and Q will be investigated. Such a decomposition al- 
lows measuring the contribution of separate scales to the 
transfer function T(k, i). The precise specification re- 
quires a few steps that are presented next. First, we 
define the truncated velocity field as 



>(M) 



u a (k,t) ifk e Porks 
otherwise. 



(16) 



Based on this truncated velocity field we may compute 
the energy transfer involving the wavevector k and all 



T PQ (k,i) 



T PP (k, t) 



if P = 



-(T PQ (k,t) -Tpp(k,t) 



(17) 



-T Q Q(k, t) 



if 



+ 



where 



TpQ(M) = tt*(M)*!f'<»(k,t). 



(18) 



The nonlinear term 'Fq ' (k, t) is defined by the convo- 
lution of the truncated fields: 

*^(k,t)=M afh £ ^(p,t)<^(q,t), (19) 

p+q=k 

where the sum is over all triads with p € P and q G Q 
such that p + q = k. 

For a statistically isotropic, homogeneous turbulence it 
is convenient to average over spherical shells in wavevec- 
tor space. In addition, in view of the considerable com- 
putational effort involved in computing all interactions 
between the very large number of scales present in the 
flow, we introduced a slight coarse-graining in terms of 
the regions P and Q as shown in Fig. [5] Specifically, 
it was found adequate to group together contributions 
from four adjacent shells. Other more coarse 'groupings' 
of wavenumbers have been considered in the literature 
with the aim of extracting the dominant interaction pro- 
cesses at a reasonable computational effort. As an exam- 
ple a 'logarithmic' grouping was adopted in |16| combin- 
ing contributions from bands with a width of 2 k . In this 
paper we will look at the interactions of four shells P at 
distance k p (cf. Fig. [2J with four shells Q at distance k q 
that contribute to the nonlinear energy transfer to shell 
§/t characterized by the wavenumber k. 

In terms of the transfer function TpQ(k, i) we may now 
define the required spectral transfer functions. The en- 
ergy transfer term Q17f) gives the exchange of energy by 
the triad (k, p, q) where the latter two wavevectors are 
specified by the sets P and Q and the triangle constraint. 
Summing over all modes k in shell §& we obtain the ex- 
act exchange of energy in the fc-th shell between k, k p 
and k„: 



T pq (k,k p ,k q ,t) = V T m (k,t) 



(20) 



We refer to T pq as the 'three-mode' transfer. The total 
energy transfer function T(k, t) can be computed directly 
from © or as sum of the contributions from l|20|) : 

T(k,t)=JZ T p (k,k p ,t), (21) 

in which the 'two-mode' transfer T p is given by: 

lp(k, kp,t) = / -L pq {k, k p , k q ,t). \^£) 

*■ — 'k„ 



The individual transfer-terms T(k,t), T p (k,k p ,t) and 
T pq (k, fc p , k q , t) give respectively more detailed character- 
istics of the energy transfer. The total transfer T(k,t) 
expresses the amount of energy transferred from (nega- 
tive) or to (positive) shell §&. All three transfer-terms 
T, T p and T pq will be used to investigate the transfer of 
energy in the sequel. 



D. Simulation details 

The numerical integration of the Navier-Stokes equa- 
tions JJl is done via a four-stage, second-order, compact- 
storage, Runge-Kutta method J34J. To fully remove the 
aliasing error we applied a method that employs two 
shifted grids and spherical truncation |35f. We con- 
sider the canonical problem of forced turbulence in a cu- 
bic box of side Lb with periodic boundary conditions. 
Direct numerical simulations are characterized by N 3 
computational points, where N is the number of grid- 
points used in each direction. A detailed description 
of the simulation setup and the validation of the nu- 
merical procedure can be found in 21] . The compo- 
nents of the wavevector k are fc Q = (2t: / Lb)n a where 
n a =0,±l,±2,...,±(JV/2-l),-iV/2 for a = 1,2,3. 
The numerical simulations are defined further by the size 
of the domain (L{,—1), the computational Reynolds num- 
ber Re and the energy injection rates to the two distinct 
bands (e Wi i, e m , 2 ). 

We will study this homogeneous turbulent flow at two 
different computational Reynolds numbers, i.e., Re — 
1061 and Re = 4243. In case of homogeneous, decay- 
ing turbulence these Reynolds numbers correspond to 
R\ = 50 or 100, in terms of the initial Taylor-Reynolds 
number [2l|. The large-scale forcing of Si has an energy 
injection rate s w> i = 0.15 that is used as reference case. 
For all simulations the fractal dimension was kept con- 
stant and equal to Df = 2.6 J27|]. The smallest length- 
scale that should be accurately resolved depends on the 
size of the box, viscous dissipation and energy injection 
rate. Usually it is required that fc max ^ > 1 |M Hy, H3 
in terms of the Kolmogorov length-scale r\ and the max- 
imal magnitude of the wavevector fc max = irN/Lb that 
enters the computations. In our simulations k max rj <J 2 
indicating that the small scales are well resolved. 

We consider time-averaged properties of the turbulent 
flow. For a function h these are defined by 



(h) t 



lim 

t— »oo t — to 



h{r)dT 



T-t 



h(r)dT, (23) 



where T is sufficiently large. We start the averaging 
at to = 5 which corresponds to about 10 eddy-turnover 
times for the simulated cases. The final time was taken 
equal to T = 30, so all results are averaged over approxi- 
mately 50 eddy-turnover times. The accuracy of this ap- 
proximation to the long-time average, measured as the 



ratio of the standard deviation and the mean signal is 
less than 5% for all investigated quantities. 

The energy spectra presented in this paper are shcll- 
and time-averaged. Moreover, we focus on compensated 
spectra E c in which we use non-dimensional Kolmogorov 
units: E c (k) = (e); 2/3 k 5 / 3 (E s (k,t)) t in terms of the 
shell-averaged spectrum E s (k,t) — J2kes E(k,t). The 
compensation of the spectrum is not strictly required to 
observe the characteristic changes in the energy distribu- 
tion, but as it gives more information about the dominant 
scales present in a flow it will be used throughout. 



III. BROAD-BAND FORCED TURBULENCE 

To investigate the energy dynamics in broad-band 
forced turbulence, first the influence of variation of the 
strength and location of the second forced band K mp 
on the global characteristics and spectrum of the flow is 
presented (fill All . Then we examine in more detail the 
effects of these variations on the energy distribution and 
transfer characteristics in the system (IIIIBfl . 



A. Energy distribution in forced turbulence 

We first concentrate on the application of the high-fc 
forcing band at different locations in spectral space. We 
apply a constant energy input rate e w _2 = 0.15 to this 
band. Simultaneously, the large-scale forcing to the first 
shell §i is e Wl i = 0.15. The computational Reynolds 
number is Re — 1061. We forced the bands 



4j,P+3 



for 



p = 5, 9, 17, 25. The parameters of these simulations with 
some of the statistics are further presented in the Ap- 
pendix (Runs 1 and 14—17 in Table HI and ITU are con- 
cerned here). 

The total kinetic energy, energy dissipation rate and 
Taylor- Reynolds number are shown in Fig.|21as a function 
of the location of the left-boundary p of the high-fc forced 
band K P:P+ 3. The first data point refers to the classical 
large-scale forcing only (Run 1). Application of broad- 
band forcing in the different bands changes the charac- 
teristics of the flow modifying primarily the amount of 
small scales. This forcing in the second band is seen to 
increase the energy dissipation in the system. The Kol- 
mogorov dissipation-scale and the Taylor-Reynolds num- 
ber decrease, suggesting that the characteristic scale at 
which dissipation plays an important role is shifted to 
smaller scales. We notice that the total energy in the 
system is only slightly affected by the introduction of 
forcing. Moving the broad-band forcing to very small 
scales implies that there is no longer a strong influence 
on the flow because the energy injected in the small scales 
appears to be also dissipated immediately. 

The compensated, shell- and time-averaged, energy 
spectrum for different locations of the forced region 
IKp.p+3 is shown in Fig. 0] We may observe that the 
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FIG. 3: Time- averaged total kinetic energy E (solid), total 
energy dissipation rate e (dotted) and Taylor-Reynolds num- 
ber R\ (dashed) for forced turbulence with different locations 
of the second band at Re — 1061. 




FIG. 5: Time-averaged total energy E (solid), total en- 
ergy dissipation rate e (dotted) and Taylor- Reynolds num- 
ber R\ (dashed) for two-band forced turbulence with varying 
strength in the second band e TO ,2 at Re — 1061. 
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FIG. 4: Compensated shell- and time-averaged energy spec- 
trum E c (k) for two-band forcing: Si and K p>p +3. Large- 
scale forcing Si (solid) , additional forcing in the second band 
K 5i 8 (dashed), Kg,i2 (dash-dotted), Ki7,2o (>), K25.28 (o) at 
Re = 1061. 



forcing causes a non-local depletion in the energy spec- 
trum for the larger scales while the tail of the spectrum 
is less affected. The pile up in the energy spectrum near 
the forcing region is characteristic of the explicit high-fc 
forcing and is suggestive of a 'blocking' or reverse cascad- 
ing. If the separation between Si and the high-fc band 
is reduced, then the interaction is stronger and a consid- 
erable depletion of the energy levels in the largest scales 
arises. This is in agreement with the large-scale forced 



FIG. 6: Compensated shell- and time-averaged energy spec- 
trum E c (k) for forced turbulence in the band Ki7,2o at dif- 
ferent strengths of forcing e W: 2 and Re — 1061. Large- 
scale forcing only (solid) , additional second band forcing with 
e w ,2 = 0.07, 0.15, 0.30, 0.45, 0.60, 0.75, 0.90 denoted as □, dot- 
ted, dashed, dash-dotted, >, 0, o, respectively. In each case 
e w .i = 0.15. 



turbulence results, where the local and non-local interac- 
tions were found to be energetically dominant while the 
distant interactions were m ainly responsible for transfer- 
ring structural information 29] . 

An effective modulation of turbulent quantities is pos- 
sible not only by a change in the range of forced modes 
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FIG. 7: Value of the maximum energy _E ma x = maxt (E c (k)) 
in the forced region for four-shell Ki7,2o (°) and eight-shell 
Ki7,24 (>) forced turbulence, measured from the compensated 
spectra at Re = 1061. 



but also via a change in the energy input rate. To investi- 
gate this we adopted an energy injection rate s w> i = 0.15 
for the large-scale forcing in §j and we vary the inten- 
sity of forcing in the second band by changing £ w p. Wc 
adopted the following values for £ w p: 0.07, 0.15, 0.30, 
0.45, 0.60, 0.75 or 0.90 and considered forcing of four or 
eight shells in Ki7 :2 o or Ki7 i2 4, respectively. The param- 
eters and characteristic quantities can be found in Table 
H|and|n]as R- uns 2 - 7 and 8 - 13. 

The total energy in the system is only slightly affected 
by the forcing-strength in the second band as shown in 
Fig [SJ An increased forcing strength introduces addi- 
tional energy into the flow at small scales that is dis- 
sipated very efficiently. This is expressed by the linear 
increase in (s)t- I n Fig. we present the compensated en- 
ergy spectrum for various strengths of the forcing e Wi 2- 
The energy in the forced region reaches higher values 
with increasing s w> 2- This may be further observed in 
terms of the energy maximum E max = max^ (E c (k)) in 
Fig. [7| Changing the strength of the broad-band forc- 
ing induces a characteristic depletion in the larger scales. 
This suggests that the additional forcing term enhances 
the nonlinear interactions, which influence various scales 
quite far away from the forced region. The energy that 
is injected at the larger scales is transferred even more 
effectively through the cascade as e W: 2 increases. This 
effect appears similar to the so-calle d sp ectral short-cut 
observed in nature and experiments |25j . In case of such 
a short-cut the energy from larger scales is diverted quite 
directly to fine scales largely by-passing the traditional 
cascading. This mechanism was explained in the case of 
flow over forest-canopies in [25j . We will investigate it in 
more detail in the next section. 



FIG. 8: Location of the maximum energy in the forced region 
for four-shell Ki7,2o (°) and eight-shell Ki7 : 24 (>) forced turbu- 
lence, measured from the compensated spectra at Re = 1061. 



While the energy-maximum -E max shifts to higher val- 
ues with increasing strength of the high-fc forcing, the 
location of the peak moves towards larger scales. This 
may be seen from the value of (kr]) max at which the 
maximum of E c (k) is attained (cf. Fig. |S}. If we inject 
the same amount of energy to Ki7 j2 4 that is two times 
wider than Ki7 j2 o the maximum of the response decreases 
(cf. Fig.EJ) while the location of the peak moves towards 
smaller scales (cf. Fig. [SJ. These numerical experiments 
are in agreement with observations for decaying turbu- 
lence in which initially additional energy is assigned to 
small scales [22|. The non-local energy cascade toward 
the small scales was found to increase remarkably during 
the initial period of decay. This is quite similar to our 
observed increase of the dissipation rate arising from an 
increased forcing of the high-fc band. 

A final quantification of the non-local effect on the 
spectrum that arises from the high-fc forcing is collected 
in Fig. 03 Here we displayed the normalized accumulated 
energy 



S E (k) 



E k < <k E c (k>) 
J2 k E c (k) 



(24) 



in the consecutive shells. As pointed out, varying the 
properties of a flow in a specified spectral region can 
change the behavior of a flow well outside this region. 
In terms of S^fc) we notice that close to 90% of the en- 
ergy is present in the first 10 shells (Fig. |5J) when only 
the large-scale forcing is applied. Influencing the flow at 
smaller scales in Kj^^o is seen to remove most of the en- 
ergy from these larger scales while there is only a slight 
impact on the dynamics of small scales. This effect be- 
comes more pronounced with increasing e Wi i- The under- 
lying changes in the energy transfer will be considered in 
more detail in Sec. IIVI 
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FIG. 9: Normalized accumulated energy SE(k) for two-band 
Si and Kir^o forced turbulence at Re = 1061. Large-scale 
forcing at e W: i — 0.15 (solid) with additional second band 
forcing at e w>2 = 0.07,0.15,0.30,0.45,0.60,0.75,0.90 is de- 
noted by the □, dotted, dashed, dash-dotted, >, o and o 
curves, respectively. 







-5 



FIG. 10: Time-averaged energy transfer T(k) for two-band 
forced turbulence. Large-scale forcing Si (solid) with addi- 
tional forcing in the Kg,8 (dashed), Kg,i2 (dash-dotted), Ki7,20 
(>), K 25 ,28 (o) band at Re = 1061. 



B. Energy transfer spectra 

The transfer of energy in turbulence can be described 
in spectral space as interactions of triads of wavevectors 
(k, p,q) that form triangles, i.e., k = p + q. Direct 
numerical simulation with large-scale forcing shows that 
non-local interactions between wavevectors combine into 
a local energy flow [Ij, y(| • By applying forcing that is 
located in a high-fc range of spectral space we perturb the 
'natural' cascading process. The associated changes in 
the transfer of energy will be investigated in more detail 
in this subsection. Specifically, we focus on the energy 
transfer and energy transport power spectra. 

In large-scale forced turbulence energy is injected into 
the first shell and removed by the transfer term. This 
gives rise to negative values for the energy transfer in 
the forced region. In the higher shells the transfer func- 
tion takes on positive values which illustrates the trans- 
fer of energy through the cascade toward higher fc. By 
invoking the broad-band forcing we influence this basic 
energy cascade. This is clearly seen in the energy trans- 
fer spectrum which develops distinctive regions where 
T(k) = (T(k,t)) t is negative. In Fig. E3 the effect of 
variations in the spectral support of the forcing is shown 
while Fig. II II characterizes changes in the transfer func- 
tion due to an increased forcing strength of the high-fc 
band. The transfer function reaches lower values between 
the low- and high-fc forcing regions compared to the large- 
scale forced case. The reverse situation appears near the 
high-fc forced band where the transfer increases with an 
increase of the forcing intensity. This is in agreement 
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FIG. 11: Time-averaged energy transfer T(k) for two-band 
§i and Ki7,2o forced turbulence for different strengths of 
forcing in the second band e w> 2 at Re — 1061. Large- 
scale forcing (solid) with additional second band forcing at 
£ m ,2 = 0.07, 0.15, 0.30, 0.45, 0.60, 0.75, 0.90 denoted by □, dot- 
ted, dashed, dash-dotted, >, o and o curves, respectively. 



with the energy spectra presented earlier, where we ob- 
served the depletion of energy between the forced regions. 
This effect can be observed more directly from spectra of 
energy transport power that will be presented next. 

The energy transport power gives the rate at which en- 
ergy is transferred from shells fc' < fc to those with k' > k: 
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FIG. 12: Time-averaged transport power spectra for broad- 
band forced turbulence in the 1X17,24 band at Re = 1061. 
Large-scale forcing (solid) with additional second band forced 
at e w>2 = 0.15,0.30,0.45,0.60,0.75,0.90 denoted as dotted, 
dashed, dash-dotted, >, o, o curves, respectively. The supple- 
mentary forced region k\ < k < k-z is denoted with dotted 
lines. 
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FIG. 13: Time-averaged transport power spectra for Reynolds 
number Re = 4243 and forcing in the Ki7,20 band. Large-scale 
forcing (solid), additional broad-band forcing in the Ki7,2o 
band with e Wt 2 — 0.30 (dashed). 



II(fc) = (II(fc, t))t in Fig. E| for forcing with various 
strengths in the K47.24 band. In case of large-scale forc- 
ing only, the transport power is positive for all k as the 
energy is transferred toward smaller scales and reaches 
zero for large k indicating the general property of the to- 
tal transfer function T(t) = 0. The application of high-fc 
forcing for k\ < k < k% changes this well-known picture. 
First, we note that the values of the transport power are 
all similar in the largest scales, where the flow is gov- 
erned by the same energy input. The transport power 
for < k < k\ becomes larger at higher e w ^- A strik- 
ing change of the behavior arises for k near and inside 
the high-fc forced region. The transport power spectrum 
even assumes negative values for fc rs fci. 

The observed behavior of the transport power in 
Fig- El is partly due to the relatively low Reynolds num- 
ber that was used. At sufficiently high Reynolds num- 
bers, the dissipation scales are much more separated from 
the high-fc forced scales. In this case a plateau of IT 
will arise at low wavenumbers: H(k,t) » e Wt \ for fc low 
enough 1]. This property is not observed at the compu- 
tational Reynolds number considered so far. 

In cases specified by Runs 18 — 19 we consider the 
flow at a four times higher computational Reynolds num- 
ber. The overall results for the energy spectra and energy 
transfer were found to be qualitatively the same as in the 
lower Reynolds number cases. However, a plateau may 
now be observed in Fig. 1131 where we present the trans- 
port power for the higher Reynolds number. In this case 
the transport power does not decrease below zero in the 
forced region. The second forcing band is well separated 
from the dissipation region and the transport power in 
this band is much larger, approaching a maximum 0.23 
that is near the energy injection rate e Wj 2. 

In this section we have looked at the effect of high-fc 
modulation of the energy cascading process that leads 
to an increased energy dissipation in small scales. This 
process is supported by an increased energy transfer to 
smaller scales via nonlocal triad interactions. The effect 
of increased energy rate by the application of broad-band 
forcing is seen in the energy transfer and transport power 
spectra. In the next section we will look more closely 
at the interactions of various scales of motion under the 
influence of broad-band forcing by considering the two- 
and three-mode transfers T p and T pq introduced in i|22|) 
and J2IJ). 



IV. TWO- AND THREE-MODE INTERACTION 
OF SCALES 



^max K 

U(k,t)= J T(k',t)dk' = - jT{k\t)dk', (25) 



where fc max = irN/Lb is the cut-off wavenumber. 
We present the time-averaged transport power spectrum 



The energy dynamics of turbulent flow is generally dis- 
cussed in terms of the transfer of kinetic energy from 
larger to smaller scales through nonlinear interactions. 
The statistical properties of turbulence are determined by 
these interactions. In the previous section we have shown 
how additional broad-band forcing of inertial range scales 
can modify the classical picture of the Kolmogorov cas- 
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FIG. 14: The normalized triad energy transfer function 
T pq (k,k p ,k q ) for k/(2n) = 14,42,82 in (a), (b), (c), respec- 
tively. The dashed lines correspond to the lower (fei) and 
upper (k%) wavenumbers used in the broad-band forcing at 
Re = 4243. The contour levels are ±1/2", n = 0, . . . , 18 that 
are the same for all three pictures. 
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FIG. 15: The normalized triad energy transfer function 
T pq (k,k p ,k q ) for k/(2%) = 46 at Re = 4243. The upper 
(bottom) rectangle presents broad-band (large-scale) forced 
turbulence respectively. The contour levels are ±1/2", n = 
0,...,24. 



cade. To investigate the observed turbulence modulation 
effects in more detail we consider the underlying two- and 
three-mode energy transfer terms in this section. This 
will clarify to some extent the changes in the various 
nonlinear interactions that give rise to the observed al- 
terations in the spectra and energy transfer. 

We start with the three-mode transfer that is averaged 
in time T pq {k,k p ,k q ) = (T pq (k, k p , k q ,t))t and split this 
term into its positive and negative parts: 

l pq {k, k p: k q ) = l pq {k, k p , k q ) ± l pq {k, k p , k q ), (zo) 

in which 

l-pq\k, kp, K q ) II i pq \k, Kp, kq) > U, 



^pq\^> ^v ™q 







otherwise, 



(27) 
with a similar definition for the negative part: 

-Lpq\ki kpi kq) it lpq\K, kp, k q ) < U, 



T (k, k p , k q ) 







otherwise. 



In terms of these contributions we examine the normal- 
ized triad energy transfer 



lpq[K : /Cp, Kq) IpAfc, Kp, Kq) 



1 min \rC) 



-L max v ) 



(29) 



where 



■*minv"<J 



min (T pq {k,k p ,k q )) 



k i, . h'i, 



T max (k) = max(T p +(fc,fc p ,/c g )) 



h'p . n-q 



Through the scaling of T and T+ with T t 



pi 



(30) 



and 



T max respectively, the normalized transfer is well suited 



12 



to characterize the overall structure of the three-mode 
transfer function, even in cases in which the order of 
magnitude of T pq varies considerably. The normalized 
energy transfer T pq (k,k p ,k q ) is plotted in Fig. ITU for 
three different wavenumbers fc/(27r) = 14,42,82, based 
on Run 19 in which R\ = 75. The three fc-values that 
are selected correspond to wavenumbers below the forced 
region (k/(2ir) = 14) or to wavenumbers that are consid- 
erably larger. Such contour maps for T pq can also be 
found in [ijj for the case of large-scale forced turbulence. 
For completeness, we also presented the results from such 
large-scale forced turbulence (Run 18) comparing these 
directly to the broad-band forced turbulence (Run 19). 
This contour map is shown in Fig.^]for fe/(27r) = 46. 

The strongest interactions are observed for modes with 
wavenumbers between the largest forced scales and the 
high-fc forced region as can be seen in Fig. I14f a). As 
in the case of large-scale forcing only we observe very 
strong interactions between Fourier modes of consider- 
ably different scales. These are located in the corners 
of the rectangular domains in the k p -k q plane. Distant 
interactions are well separated from the origin in these 
figures. Their contribution to the transfer is seen to be 
very small, as also noticed earlier in the literature [19j. 
The change of sign in the transfer function that occurs at 
k p = k and k q = k respectively on the k p -k q planes indi- 
cates that in this region the energy is mainly transferred 
to higher k. 

The most efficient transfer takes place between two 
wavevectors of similar size and one of quite different size 
as seen in the corners of the rectangular area in Fig. 1141 
This is in agreement with previous numerical experiments 
reported by various authors Q,[l£}[lj|. However, com- 
pared to the case of large-scale forcing only, we now ob- 
serve quite extended, highly energetic interactions with 
the high-fc forced region. The second forced band causes 
regions with high intensity of interactions to be much 
wider compared to the case of large-scale forcing only. 
This is visible directly in Fig. 1151 The regions with posi- 
tive and negative transfer are extended from the corners 
to the wavenumber regions where the actual application 
of forcing in the second band occurs. The energy is ex- 
changed predominantly between scales that are more sep- 
arated than in case of the large-scale forced flow where 
the dominant interactions occur only in the corners. This 
is a clear indication of the stronger non-local interactions, 
mentioned earlier. 

For further clarification of energy transfer processes 
wc turn to the time-averaged two-mode energy transfer 
T p (k,k p ) = (T p (k,k p ))-t, which gives information about 
the interactions involving a sum over all k q wavenumbers 
at fixed k and k p . The sum involves all k q wavenumbers 
that are constrained by the triadic interactions, i.e., their 
length may vary between |k — p| and |k + p|. We normal- 
ized the two-mode transfer function T p (k, k p ) in a similar 




40 60 

FIG. 16: The normalized triad energy transfer function 
T p (k,k p ) at Re = 4243 (Run 19). The contour levels are 
±1/2™, n — 0, . . . , 18. The supplementary forced region 
fei < k < &2 denoted with dashed lines. 




FIG. 17: The triad energy transfer function T p (k,k p ) at 
k p /(2n) = 30, 34, . . . , 94 and Re = 4243 (Run 19). 



manner as T p (k,k p ,k q ): 



m ( h , \ _ T p (fc, k p ) T p (fc, kp) 
±p{K,Kp) - + , 

J- mm V / maxV' v y 



(31) 



where Tp, T min and T max are defined in terms of T p 
in a manner analogous to the definitions in (|26|) . H27fl . 
(|28|l and (|30[) . In Fig. Elwe plotted the contour map 
of T p (k,kp). For larger wavenumbers this quantity was 
found to look quite similar to the case of large-scale forced 
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FIG. 18: The triad energy transfer function T p (k,k p ) 
at k p /(2ii) = 30 and Re — 1061. Large-scale forc- 
ing (solid) with additional second band forced at e w ,2 = 
0.07,0.15,0.30,0.45,0.60,0.75,0.90 denoted as □, dotted, 
dashed, dash-dotted, >, o, o, respectively. 



turbulence. The two-mode transfer function changes 
sign from negative to positive at k — k p indicating a 
downward energy flow. Comparing this to the large- 
scale forced turbulence we observe (i) strong influence 
of forcing in the regions where it is applied (denoted 
with dashed lines), (ii) extended negative energy trans- 
fer region with comparatively high magnitude above the 
k = k p line, (Hi) amplification of the backward energy 
transfer indicated by the positive region for small k and 
large k p . This region is separated from the intense nega- 
tive energy transfer region by the indicated accumulation 
of contour lines above the k = k p line appearing as the 
curved black line. 

A more quantitative overview is plotted in Fig. 1171 
displaying the two-mode transfer function in the range 
k p /(2ir) = 30, 34, . . . , 94. This clearly shows the cascad- 
ing character of the energy flow from larger to smaller 
scales in the system. The modification due to the high- 
k forcing expresses itself by the sequence of one slightly 
positive, two quite negative and one quite positive lo- 
cal extrema. The intensity of the energy transfer de- 
creases with increasing wavenumbers as less energy needs 
to be transferred. This corresponds directly to the mag- 
nitude of T m i n (fc) and T max (fc) used in the normalization 
of T p (k, k p ) (|31|) . The part in which the transfer is neg- 
ative is much wider in the broad-band forced case com- 
pared to the large-scale forced turbulence results. 

We conclude by considering the effect of varying the 
forcing strength e Wt 2 at a characteristic wavenumber 
k p /(2ir) = 30 on the two- mode energy transfer function 



T p (k,kp). This is shown in Fig. ^J In the large-scale 
forced case at e Wi % = the transfer is very small com- 
pared to the cases in which the high-fc forcing is active. In 
addition, the effect is very localized (solid line in Fig. lTH)) . 
The forcing in the high-fc band completely changes this 
behavior. The intensity of the energy transfer is directly 
related to the value of s w> 2- Additional extrema appear in 
the two-mode transfer function. The high-fc forced cases 
display two pairs in which a negative minimum is com- 
bined with a positive maximum, while large scale forcing 
only yields one such combination. Correspondingly, the 
min-max pair at high fc is associated with the large scale 
forcing in Si while the min-max pair at lower k origi- 
nates from the additional forcing in the second band. We 
also investigated three-band forcing and observed further 
peaks in the energy transfer spectra. 



V. CONCLUDING REMARKS 

We performed direct numerical simulations of broad- 
band forced turbulence to explore accumulated effects on 
the time-averaged energy transfer in isotropic homoge- 
neous turbulence. Using broad-band forcing based on a 
recently proposed mathematical model for a fractal stir- 
rer 123 we have shown how the application of such forc- 
ing modulates turbulence both qualitatively and quan- 
titatively. The modulation is similar to that observed 
in experiments based on flows through porous media or 
canopies. Specifically the perturbation of a flow arising 
from the contact with complex physical boundaries en- 
hances the dissipation and causes an abrupt energy drain 
from large to small scales. This aspect of simultaneous 
perturbation of a flow on a spectrum of length-scales is 
retained in the cases studied here. 

We found that broad-band forcing that perturbs a tur- 
bulent flow at smaller scales enhances non-local triad in- 
teractions and alters the detailed cancellation processes 
that occur in the traditional large-scale forced flows. This 
leads to non-local modifications in the energy transfer 
spectrum and the energy distribution among scales. We 
verified this by partitioning the nonlinear term in the 
Navier-Stokes equations in terms of different triad contri- 
butions to the total transfer function. The energy trans- 
port power is found to be enhanced in the spectral region 
in between the large-scale and the high-fc forced bands. 
This characteristic may be influenced via the control pa- 
rameters of the applied forcing, i.e., its strength and ex- 
tent of agitated scales, and allows optimizing transport 
processes of turbulent flows. 

Future study will involve the examination of the con- 
sequences of forcing in the physical space context. We 
will investigate the geometrical statistics of broad-band 
forced turbulence looking at the interactions of strain and 
vorticity and their modulation by the applied forcing. 
This may help understanding which physical processes 
are responsible for the observed modulations and how to 
exploit this to enhance physical space mixing. 



14 



Run 


&w 


m 


P 


Run 


£w 


771 


P 


1 


0.15* 


- 


- 










2 
3 

4 
5 
6 

7 


0.07 
0.30 
0.45 
0.60 
0.75 
0.90 


17 
17 
17 
17 
17 
17 


20 
20 
20 
20 
20 
20 


8 
9 
10 
11 
12 
13 


0.15 
0.30 
0.45 
0.60 
0.75 
0.90 


17 
17 
17 
17 
17 
17 


24 
24 
24 
24 
24 
24 


14 
15 
16 
17 


0.15 
0.15 
0.15 
0.15 


5 
9 

17 
25 


8 
12 

20 

28 


18 

19 


0.15* 
0.30 


17 


20 



TABLE I: Direct numerical simulation parameters using a res- 
olution of N = 128 and Re = 1061 in Runs 1 — 17, and a res- 
olution of N = 256 at Re = 4243 in Runs 18 - 19. The cases 
with large-scale forcing only are denoted by *. In this table 
e w denotes the energy input-rate in the high-fc band, except 
Run 1 and 18 in which it corresponds to the energy input-rate 
in Si. Moreover, m and p characterize the spectral support 
of the high-fc band K m , p . 



APPENDIX 

The main parameters of the simulations are collected 
in Table [I] The corresponding statistics of the veloc- 
ity fields are summarized in Table [n] The quanti- 
ties compiled in Table [ij are the Kolmogorov dissipa- 
tion wavenumber kd which is the inverse of the Kol- 
mogorov length-scale r\, the product fc max ?7, the Tay- 
lor microscale A = (5E/J2 k k 2 E(k, t)) 1 / 2 , the Taylor- 
microscale Reynolds number R\ = Xu'/v, the inte- 
gral length-scale L = 3vr/(4E) £ k fe -1 .E(k, t), the in- 
tegral Reynolds number R L — Lu'/v, the r.m.s veloc- 
ity u' = (2E/3) 1 / 2 , the energy dissipation rate s = 
^ k 2^fc 2 _E(k, t), the eddy-turnover time r = L/u' and 

the skewness S = 2/35 (A/V) 3 £ k k 2 T{k,t). All these 
quantities in Table ITU are time- averaged (•)( as described 
in Sec. IlTdI 

We also checked that the alteration of the cascading 
process caused by the high-fc forcing does not influence 
the isotropy of the flow field. A measure of isotropy was 
suggested in [38| given by: I 2 (t) = ij)\(t) / 'ijJiif) where 
M*) = <| ei (k)u(M)| 2 ), Mk,t) = <|e 2 (k)u(k,i)| 2 ) are 



the kinetic energy along the components of two orthogo- 
nal solenoidal unit vectors ei(k) = k x z(k)/|k x z(k)|, 
e2(k) =kx ei(k)/|k x ei(k)| where z(k) is a randomly 
oriented unit vector. The operator (■) denotes averaging 
over these random unit vectors. For isotropic turbulence 
one can expect to find 1 = 1, i.e., ip\ = -02 which was 
confirmed to close approximation in all simulations. De- 
viations from the expected value for / were found to be 
of the order of 1%. 



Run 


kd 


"'max'] 


A 


Rx 


L 


Rl 


1 
u 


e 


T 


S 


1 


116 


3.26 


0.123 


52 


0.23 


97 


0.40 


0.15 


0.57 


0.49 


2 
3 

4 
5 
6 

7 


130 
156 
168 
178 
186 
193 


2.91 
2.42 
2.25 
2.12 
2.03 
1.95 


0.100 
0.069 
0.061 
0.056 
0.051 
0.048 


43 
30 

27 
25 
23 
22 


0.23 
0.23 
0.23 
0.22 

0.22 
0.22 


100 
98 
101 
102 
102 
100 


0.41 
0.41 
0.42 
0.43 
0.43 
0.43 


0.24 
0.50 
0.66 
0.83 
1.00 
1.16 


0.56 
0.54 
0.54 
0.52 
0.52 
0.50 


0.35 
0.21 
0.17 
0.15 
0.14 
0.13 


8 
9 

10 
11 
12 
13 


140 
156 
169 

178 
187 
194 


2.68 
2.41 
2.24 
2.11 
2.02 
1.94 


0.085 
0.070 
0.060 
0.054 
0.049 
0.046 


37 

31 
26 
24 
22 
21 


0.23 
0.23 
0.22 

0.22 
0.22 
0.21 


99 
102 
98 
100 
97 
96 


0.41 
0.41 
0.41 
0.42 
0.42 
0.42 


0.33 
0.50 
0.68 
0.85 
1.03 
1.20 


0.56 
0.56 
0.54 
0.53 
0.52 
0.50 


0.29 
0.22 
0.18 
0.16 
0.15 
0.13 


14 
15 
16 
17 


138 
138 
140 
143 


2.73 
2.72 
2.69 
2.64 


0.090 
0.089 
0.084 
0.082 


39 
39 
36 
36 


0.22 
0.23 
0.22 
0.23 


96 
102 
96 
98 


0.41 
0.42 
0.40 
0.41 


0.30 
0.31 
0.33 
0.35 


0.52 
0.56 
0.55 
0.56 


0.44 
0.39 
0.28 
0.22 


18 
19 


325 
432 


2.32 
1.74 


0.065 
0.040 


115 
75 


0.21 
0.21 


368 
394 


0.42 
0.44 


0.15 
0.46 


0.50 
0.47 


0.51 
0.38 



TABLE II: Direct numerical simulations statistics of the dif- 
ferent cases studied. 
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